Coherent electron transport by adiabatic passage in an imperfect donor chain 
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Coherent Tunneling Adiabatic Passage (CTAP) has been proposed as a long-range physical qubit 
transport mechanism in solid-state quantum computing architectures. Although the mechanism can 
be implemented in either a chain of quantum dots or donors, a ID chain of donors in Si is of particular 
interest due to the natural confining potential of donors that can in principle help reduce the gate 
densities in solid-state quantum computing architectures. Using detailed atomistic modeling, we 
investigate CTAP in a more realistic triple donor system in the presence of inevitable fabrication 
imperfections. In particular, we investigate how an adiabatic pathway for CTAP is affected by 
donor misplacements, and propose schemes to correct for such errors. We also investigate the 
sensitivity of the adiabatic path to gate voltage fluctuations. The tight-binding based atomistic 
treatment of straggle used here may benefit understanding of other donor nanostructures, such as 
donor-based charge and spin qubit s. Finally, we derive an effective 3x3 model of CTAP that 
accurately resembles the voltage tuned lowest energy states of the multi-million atom tight-binding 
simulations, and provides a translation between intensive atomistic Hamiltonians and simplified 
effective Hamiltonians while retaining the relevant atomic-scale information. This method can help 
characterize multi-donor experimental structures quickly and accurately even in the presence of 
imperfections, overcoming some of the numeric intractabilities of finding optimal eigenstates for 
non-ideal donor placements. 

PACS numbers: 05.60.Gg, 73.63.Kv, 73.21. La 



I. INTRODUCTION 

Silicon quantum bits (qubits) are pursued due to the 
promise of long spin coherence times and the processing 
expertise of the semiconductor industry, that potentially 
could produce a quantum computer (QC). Over the past 
decade, many different Si qubits have been proposed. 
Some of these are based on nuclear [1] or electronic spins 
0,0] of donors, while others are based on electronic spins 
in lithographically or electrostatically defined quantum 
dots 0, and still others are based on the localized 
charge states @ of confined electrons. 

QC architectures will benefit greatly if long-range co- 
herent quantum transport schemes can be incorporated 
into the architecture for efficient transfer of the qubit 
state to different areas of the QC. In Ref 0], a 2D bilin- 
ear array architecture for donor qubits was presented that 
utilized a novel non-local quantum transport mechanism 
called coherent tunneling adiabatic passage (CTAP) jgj]. 

In essence, CTAP is a solid-state analogue of the well 
known STIRAP (stimulated raman adiabatic passage) 
protocol of quantum optics [9j. It is ideally suited to 
physically transporting quantum information across a 
chain of donors or quantum dots. STIRAP has been re- 
cently demonstrated in optical waveguides using photons 
I3-E3I- A P art from dot SI and donor [7] based CTAP 
for quantum information transport, CTAP has also been 
proposed for transporting single atoms [l3|, [li| and Bose- 
Einstein condensates 15, [l6[. In previous works, com- 



parisons have been made between the solid-state and 
quantum optics versions of CTAP 17, 18[. Recent works 
have explored alternate coupling schemes [HI and CTAP 
based interferometry [2o[ in donor nanostructures. Since 
triple dot structures have already been fabricated 
[25[, solid-state CTAP may not be far from reality. 

In a system of tunnel-coupled donors or quantum dots, 
CTAP can be realized by modulating barrier gates (Fig. 
[TTa)) adiabatically in a counter intuitive pulsing sequence 
[8[ . The protocol begins with an electron localized at one 
end of the chain in a superposition of its up and down 
spin states, and ends with the electron transported to the 
other end of the chain while retaining the quantum super- 
position in the spin basis (Fig. Gib)). At all times during 
the transfer, there is negligible occupation probability of 
the electron at any point along the chain except for the 
ends. In effect, the protocol realizes certain pathways in 
the eigenspace through sensitive control of the molecular 
states of the donors or the dots. The method is expected 
to reduce gate densities relative to other quantum trans- 
port architectures in QCs, such as electron shuttling on 
the surface [261 ] . 

Bottom-up approaches have been successful in placing 
single donors in Si with a precision of 1 nm [27j . Recently, 
a ID wire of donors has been fabricated by STM pat- 
terning single Phosphorus donors in Si [28]. Other recent 
STM patterned structures include a 2D delta doped layer 
of P donors for gating other nanostructures [29] , a single 
P donor surrounded by four leads [3pj , and a quantum 
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FIG. 1: (a) A three-donor CTAP device with surface (S) gates 
and barrier (B) gates. The system has 3 ionized donors (L 
= Left, M = Middle, R=Right) and 1 bound electron. The 
device size is 60 nm x 30.4 nm x 30.4 nm. (b) The localization 
of the donor electron at 3 different voltage configurations in 
their stepped sequence. Although CTAP is a charge transport 
mechanism, in QC architectures it can be used to transport 
the spin superposition from the left to the right donor. 

dot formed by patterning a few donors [3l|. Top-down 
approaches by selective ion-implantation have also been 
successful in fabricating few donor devices (32|. In Kef 
[33( | , strategies to build a few donor CTAP device by ion- 
implantation was explored. Such experimental progress 
in donor placement motivates the examination of donor- 
CTAP's robustness to non- ideal placement of the donors 
and voltage fluctuations. 

Since the inter-donor tunnel coupling is sensitive to 
the relative locations of the donors, even small place- 
ment errors of a lattice constant or less are likely to af- 
fect the adiabatic pathway in practical implementations 
of CTAP. In fact, this donor straggle problem occurs in 
all donor qubits, whether spin or charge based. In the 
Kane qubit [l[ , the inter-donor exchange energy has been 
shown to oscillate with relative donor separations [34| . In 
a donor charge qubit, the symmetric- anti-symmetric gap 
(As as) a l so exhibits oscillations as a function of donor 
separation and orientation [35j. The sensitivity of these 
parameters occur both due to the exponential fall-off of 
the envelope wavefunction and due to the rapidly oscillat- 
ing Bloch functions resulting from the momentum states 
near the conduction band valleys of the host. A reso- 
lution to this problem is to experimentally characterize 
each donor qubit and their interactions [36-38], and ob- 
tain voltage pulses that correct for these straggle effects. 

In a previous work [39j, the existence of an adiabatic 
pathway in an ideal triple donor chain was established 
using large scale atomistic tight-binding simulations. For 



the specific donor configuration in consideration, a set 
of voltages was found to implement a complete transfer 
sequence of the electron from one end of the chain to the 
other, thereby proving that a translation of this protocol 
from the quantum optics framework to a real solid-state 
nanostructure is valid and feasible. 

In this work, we investigate the robustness of the adi- 
abatic path against imperfections such as donor straggle 
and gate voltage fluctuations. We also develop a general 
model that provides improved insight about the tunnel 
couplings and on-site energies of the multi-million atom 
system. The model can assist experiments to characterize 
donor nanostructures efficiently, and can help to design 
voltage pulses to correct for straggle. Although we re- 
strict our attention to a triple donor chain, the methods 
and the results presented here can be easily translated 
to other 3D confined nanostructures in silicon and longer 
CTAP chains. 

We have employed atomistic tight-binding (TB) the- 
ory [40|| , as this technique provides a highly accurate de- 
scription of impurities in silicon 4ll444l and allows the 
fast solution of million atom systems [45[. This method 
also treats the full band structure of the host, and pro- 
vides a unified framework to treat realistic geometries, 
gate voltages, and disorder in an atomistic setting [46 ]. 

This paper is organized as follows. In Section II, we de- 
scribe the specifics of the nanostructure considered here, 
and introduce a solution for the adiabatic pathway with 
ideal donor placements [39| . Section III elaborates on the 
method of calculation. Section IV describes the results 
in detail. We show direct TB simulations of the effect 
of straggle and its corrections. We also show an effec- 
tive 3x3 matrix model based on TB wavefunctions that 
can capture the effect of the gates and the details of the 
straggle. 



II. ADIABATIC PATH IN AN IDEAL TRIPLE 
DONOR CTAP DEVICE 

The CTAP device used in this work is shown in Fig. 
[D[a). For consistency, we have used the same device as in 
Ref [39|, except displacing the donor positions to study 
effects of straggle. The structure consists of about 3.5 
million Si atoms with one electron bound across three 
ionized P donors. Two symmetry gates and two barrier 
gates, each of 10 nm width, are placed above a 5 nm 
thick oxide layer on top of the Si lattice. Ideally, the 
barrier (B) gates modulate the tunnel barriers between 
the donors, whereas the symmetry (S) gates detune the 
energies of the end donors. The donors are buried 15 nm 
below the oxide, and are also placed 15 nm apart from 
each other in the [100] direction in the ideal case. The 
closeness of the gates in this particular device gives rise 
to significant crosstalk. Further details of the structure 
can be found in Ref 1391. 
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The three donor system can be conceptually under- 
stood with an effective 3x3 Hamiltonian, 



t LM 

tLM tMR 
t M R 



(1) 



where tLM and tMR are the gate voltage tuned tunnel 
couplings between the donors L, M, and R. We as- 
sume a simplified case for this discussion such that the 
donor on-site energies are aligned, which can be practi- 
cally achieved via tuning by the symmetry gates Si and 
S 2 . The lowest states of this Hamiltonian are, 



\* j )=a j \L)+p j \M)+ lj \R) 



(2) 



where <x/, and jj are the coefficients of the donor 
states |L), |M), and \R), respectively for the eigenstate 
j. Here, j runs from 1 to 3, with j = 1 representing the 
ground state, a, /3, and 7 are each functions of tLM and 
tMR- For ideal CTAP f} 2 = at all times. Notionally 
CTAP operates through maintaining the charge degree 
of freedom in the first excited state, ^2), and adiabati- 
cally evolving ^2), such that the wavefunction is trans- 
ferred from \L) at the start to \R) at the end through 
voltage tuning such that 72 = at the start and a 2 = 
at the end. The first excited state ^2) is protected from 
^1) and ^3) by energy gaps A12 and A23, respectively. 
We will discuss later in this paper a method to project 
the multi-million atom Hamiltonian on to a 3 x 3 sub- 
space, which can provide an accurate description of the 
full Hamiltonian for the relevant states. This model is 
one of the important contributions of this paper. 

In Fig. [2j we show the molecular states of the triple 
donor system at the midpoint of the adiabatic transfer 
path for the CTAP protocol. El represents the ground 
state, while E2 and E3 are the first and second excited 
states, respectively. CTAP transport takes effect through 
the state E2. In an ideal scenario, the middle donor be- 
tween the barrier gates can be replaced by a chain con- 
taining an odd number of donors, and the same transfer 
protocol can still hold, without the need for any addi- 
tional gates [Hj. 

In this work we will use the relative electron density in 
the center donor in the state E2 (^2)) as a metric of the 
efficacy of the CTAP transfer. The justification for this 
metric directly comes from earlier works 

@,[TMl8|. Adi- 
abatic transfer is realized through the wavefunction sym- 
metries described earlier in this section and also shown 
in Fig. [2] for the midpoint of the adiabatic path. The 
vanishing center donor density in the state E2 (^2)) is 
a defining feature of this adiabatic transfer. Any marked 
build up of center donor density in E2 indicates a devia- 
tion of the wavefunction symmetries required for CTAP, 
and is an indication that the transfer is becoming more 
non-adiabatic. 
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FIG. 2: (a) The lowest three eigenstates of the system near 
the mid-point of an adiabatic path for CTAP. These were 
obtained from tight-binding simulations of the device con- 
taining about 3.5 million atoms. The transfer takes effect 
through the first excited state of the system (E2), in which the 
electron density diminishes at the center donor. This point 
in the adiabatic path was realized for the gate configuration 
(Vsi, Vbi,Vb2, Vs2) = (-0.01325, -0.11, -0.11, 0) V. 
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FIG. 3: Population evolution through CTAP sequence, (a) 
For these calculations we use a Gaussian pulse sequence. For 
the intuitive direction, tLM is the leading pulse, followed 
by tMR- In the counter-intuitive direction, this order is re- 
versed, (b) Final state population at donor R (prr) for 
increasing total time of the pulse sequence in the counter- 
intuitive direction. Note the smooth, monotonic improvement 
in transfer fidelity with increasing time, (c) Population as a 
function of time through the CTAP protocol for total time 
T max = 207r/Q m ax- Observe the smooth population transfer 
with negligable population at the central donor, (d) As (b) 
but for the intuitive pulse sequence. Note the oscillating pop- 
ulation, with the final state population depending sensitively 
on the total time, (e) As (c) for the intuitive direction, again 
showing that the population exhibits nonadiabatic oscillations 
through the transfer protocol. 



In Fig. [3j we performed a time dependent analysis 
based on the effective 3x3 CTAP density matrix to 
show the relation between the donor densities in the adi- 
abatic and non-adiabatic regimes. The time evolution of 
the densities was obtained from dp/dt = — (i/K)[H,p\. 
We ignored effects of decoherence for this analysis as it 
is beyond the scope of this work, and was considered in 
Ref 471 in the more general context of adiabatic infor- 
mation transport. To effect an electron transfer from 
donor L to donor R, CTAP relies on a counter-intuitive 
pulsing (CIP) of the barrier gates, in which B2 is pulsed 
before B\. On the other hand, the more conventional 
sequential transport from donor R to M and then to L 
is non-adiabatic in nature, and can be realized with an 
intuitive pulsing (IP) scheme, in which B\ is pulsed be- 
fore B2. Fig. [3] shows the Gaussian pulses (in (a)) and 
the evolution of the donor densities in time (in (b), (c), 
(d), and (e)). The IP sequence produces an oscillating 
donor density in the end donor R (Fig. EJd)) and sig- 
nificant population in M (Fig. Efe)) which oscillates in 
time. This pulse does not guarantee a robust final trans- 
fer to the end donor. The CIP scheme evolves the system 
adiabatically as the density in R (Fig. [3fb)) smoothly in- 
creases to 1 with vanishing density at M (Fig. SJc)) at 
all times, making CTAP a robust transport method. 

The build up of center donor density in E2 represents a 
deviation from ideal CTAP. Imperfections such as deco- 
herence, straggle, and voltage noise can enhance this cen- 
ter density, taking the system into a more non-adiabatic 
transfer regime. The center donor density is thus a de- 
gree of measure of the non-ideality of a CTAP transfer, 
as suggested also in Refs 



III. METHOD 

Since in principle CTAP is an adiabatic problem, we 
analyze it by solving the time independent Schrodinger 
equation at different bias points. We have employed 
the 10 band sp 3 d 5 s* tight-binding model with nearest- 
neighbour interactions. The model parameters were op- 
timized by a genetic algorithm with appropriate con- 
straints to reproduce the important features of the bulk 
bandstructure of the host [40[. The model parameters 



have been well-established in the literature [45|, |49| , and 
calculations performed with these parameters have been 
verified against experimental measurements in a number 
of works gjfij m Hot 

Each P donor was modeled by a Coulomb potential 
screened by the dielectric constant of Si. At the donor 
site, a cut-off potential Uo was used, and its value opti- 
mized so that the ground state binding energy of -45.6 
meV was obtained for a donor in bulk Si. In this TB 
model, the valley-orbit interaction that lifts the six-fold 
degeneracy of the donor ground state is inherently in- 
cluded jHij. 



The electrostatic gate potential was obtained from a 
commercial Poisson solver [52] for a single gate, and the 
potential for the three gates was treated as a superpo- 
sition of the single gate solution. Although the voltages 
presented here will have some offsets from the realistic 
case due to this, the basic trends and analysis presented 
here are general, and this approach reduced the compu- 
tational resources necessary to complete the calculations 
in a tractable time. 

The net potential was then interpolated onto the 
atomistic grid for the tight-binding simulations. Closed 
boundary conditions with a model of dangling bond pas- 
sivation was used to model the interfaces [53]. The full 
Hamiltonian of about 3.5 million atoms including the four 
gate potentials was solved by parallel Lanczos and block 
Lanczos algorithms to capture the relevant eigenvalues 
and wave functions. Typical computation time for a sin- 
gle time-step of the protocol over 6 states was 7 hours on 
40 processors [54j . 



IV. RESULTS AND DISCUSSIONS 

A. Binding energy of P 3 2+ with straggle under zero 
gate bias 

In a single P donor, the lowest manifold consists of 
6 Is type orbital states, arising from the six- fold degen- 
erate conduction band valleys of Si. Due to the valley- 
orbit interaction caused by central cell effects [55], the 
six Is states are split into an orbital singlet A±, an or- 
bital triplet T2, and an orbital doublet E [56]. The A\ 
state for a P donor is separated by 11.7 meV from the 
T2 states. In a triple donor chain with donor separations 
much larger than the donor Bohr radii, the El, E2, and 
E3 states arise from linear combinations of the A± states 
of each impurity. The 12 meV gap between the A\ and 
the T2 states of a single donor also gives rise to a sig- 
nificant energy gap between the lowest 3 states and the 
higher manifold of the molecule, which is a desir- 
able condition for efficient CTAP transfer. Translation 
of this concept to quantum dots will be straightforward 
only if the valley splitting is large such that the lowest 
three states involved in CTAP are sufficiently isolated in 
energy from the valley split manifold. 

In Fig. 21 we show the lowest energy states of the P| + 
molecule at zero gate bias as a function of non-ideal place- 
ment of the middle donor. AR represents the distance 
the middle donor is displaced towards the right donor rel- 
ative to its ideal (equidistant) location. For example, in 
Fig. Ufa), AR = 0.543 nm means that the middle donor 
is displaced by one lattice constant ao (0.543 nm) along 
the donor chain. 

The states E2 and E3, which are degenerate at AR = 
0, split due to asymmetric tunnel coupling to the left and 
right donors, as shown in Fig. ffl^a), 4(b) and 4(c). The 
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FIG. 4: The lowest few states of the triple donor system at 
zero gate bias. The binding energy is expressed relative to 
the conduction band minima of Si. AR is the distance the 
middle donor is moved towards the right donor from the ideal 
(equidistant) case of 15.2 nm separation. The states for AR 
along (a) [100] (x), (c) [110] (x=y), (e) [010] (y), and (f) [111] 
(x=y=z). Plots (b) and (d) compare the lowest two (bonding 
and anti-bonding) states of a donor charge qubit (P^~) for 
donor separations along [100] (x) and [110] (x=y) respectively. 

order of meV splitting for the range of donor separations 
considered here is comparable to the order of gate voltage 
modulation for electron transfer. Hence, donor straggle 
is likely to affect the adiabatic transfer. For large donor 
misplacements of 10ao or more, there is the added prob- 
lem that the higher manifold of states move toward the 
2nd excited state E3, an undesirable condition for effi- 
cient transfer. This represents a potentially significant 
challenge to implement buried CTAP. 

Fig. 2](c) and 4(f) are for donor straggle along [110] and 
[111] directions respectively relative to the ideal chain in 
[100]. The striking feature is that the variations are al- 
most smooth, compared to the oscillations in symmetric- 
antisymmetric energy states (As as) for a charge qubit of 
P<2~ molecule shown in Fig. |H(d) and also in Ref [35j. In 
the triple donor system, the oscillations in energy states 
with donor separations seem to be somewhat mitigated. 
This indicates that relatively smooth voltage corrections 
are possible. 

Fig. 2Je) shows the effect of straggle on the donor 



spectrum for a displacement perpendicular to the chain 
along [010] direction. In this case, the spectrum is less 
affected as the middle donor remains equidistant from 
the left and right donors. The states E2 and E3 remain 
degenerate throughout as expected. 
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FIG. 5: Plots of the 1st excited state (E2) electron density for 
donor straggle of (a) AR — lao[100](x) (ao = 0.543 nm being 
the lattice constant of Si), (b) AR = 2a [100](x), (c) AR = 
la [010](y), (d) AR = 2a [010](y), (e) AR = la [100](x) 
with a corrective bias of (-0.01325, -0.11, -0.11, -0.0084) 
V, and (f) AR = 2ao[100](x) with a corrective bias of 
(-0.01325,-0.11,-0.11,-0.0170) V. Cases (e) and (f) are 
the same as cases (a) and (b) respectively, except a corrective 
voltage is applied to the S2 gate. 



B. Effect of straggle in donor position on the 
adiabatic path 

To investigate the effect of straggle more exactly, we 
have chosen the midpoint of the adiabatic path shown in 
Fig. [2] for the ideal case. Under the same gate bias, we 
have displaced the middle donor to neighboring lattice 
sites, and investigated the effect on the CTAP states. In 
particular, we compared the 1st excited state E2 obtained 
in this manner with the ideal E2 shown in Fig. O 

In Fig. 03(a) and 5(b), the middle donor is displaced by 
lag and 2ao respectively towards the right donor (in [100] 
(x) direction). It is observed that the population becomes 
dominant at the center donor. If the donor is displaced 
perpendicular to the chain ([010] (y) direction) as in 5(c) 
and 5(d), the center donor density is still negligible. This 
suggests an intuition that straggle effects are maximum 
when donors are misplaced along the chain, since the 
tunnel coupling is asymmetrically affected the most in 
such cases. 

In Fig. [5](e) and 5(d), we have found voltages that 
can be applied to the £2 gate to correct for the lao and 
2ao straggle effects of 5(a) and 5(b), respectively, even 
in the presence of considerable perturbation from other 
neighboring gates. This means that external macroscopic 
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voltages can be used to compensate for straggle effects 
due to single atom placement errors, which is a central 
point of this paper. 
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FIG. 6: Plots of the 1st excited state (E2) electron density for 
donor straggle of (a) AR = 0.5a [110], (b) AR = la [110], 
(c) AR = -0.25a [lll], and (d) AR = 0.75a [lll]. 

Similarly, in Fig. O we investigated the effect of donor 
displacement along [110] in (a) and (b), and along [111] 
in (c) and (d). These directions are expected to be more 
susceptible to rapid spatial oscillations of the host Bloch 
functions, as seen in two-donor charge qubits [35|. How- 
ever, it is found that the center population is still some- 
what negligible, which suggests that the effect of the crys- 
tal momentum states is not significant in this triple donor 
system. 

C. Sensitivity of the protocol to gate voltage 
fluctuations 



from charge traps or from small fluctuations in the gate 
pulses. Such variations are likely to affect both the tunnel 
barriers and the detuned energies of the donors, thereby 
introducing a finite population in the center donor. 

In Fig. [7J we have again used the center donor popula- 
tion as a measure of the effectiveness of the transfer, and 
investigated the robustness of the adiabatic path to volt- 
age fluctuations. For simplicity, we have modeled the 
voltage noise by perturbing Vbi by a percentage of its 
value at the mid point of the adiabatic path. 

Although Fig. [71 shows that such voltage offsets indeed 
introduce some population at the center, it is to be noted 
that center donor density remains considerably smaller 
for the more negative bias range. At -5 percent change 
in Vbi, the population at the center is still about 0.04 out 
of 1. However, the center donor density is very sensitive 
to the more positive bias range. The transfer efficiency 
is seriously hampered after +2-3 percent fluctuation in 
Vbi- 

This asymmetry in sensitivity occurs due to the order- 
ing of the CTAP states. As shown in Fig. [2j the state 
with the most central donor density is the 2nd excited 
state (E3) of the system. However, at zero gate bias, this 
state is the ground state of P| + due to strong tunnel cou- 
pling between the donors. Negative barrier gate biases 
were required to rearrange the states to realize the adia- 
batic path. Hence, more positive gate biases causes this 
state to be either the 1st excited state or ground state 
of the system, and thus breaks down the adiabatic path. 
Negative bias fluctuations, however, keep the ordering of 
the states intact, and have less effect on the path. 
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FIG. 7: Sensitivity of the adiabatic pathway to fluctuations 
in the barrier gate voltage Vbi- The plot is for the mid point 
of the path with (-0.01325, -0.11, -0.11, 0) V for which the 
center donor population is 5.4e — 4 out of 1. The center donor 
density is plotted as a function of percentage variation of \Vbi \ 
for a range of -5 to +5 percent for the ideal case as well as 
straggle of 1 and 2 lattice constants along the chain. 

The transfer efficiency of the CTAP protocol in realis- 
tic systems may be susceptible to voltage noise. In par- 
ticular, voltage fluctuations in the system can originate 



D. Effective 3x3 model developed from 
tight-binding 

Problems such as CTAP and STIRAP are well- 
described by small scale effective models, portrayed as 
a 3 x 3 matrix in the 3-donor CTAP, for example. How- 
ever, a complete and realistic description of the system 
has to include atomic scale effects from a method such 
as tight-binding. In this section, we describe a procedure 
to translate a complex atomistic system into a simpler 3 
x 3 matrix model that retains the relevant atomic scale 
information. Such a model is general beyond CTAP, and 
helps to bridge the gap between intensive atomistic the- 
ories and reduced order envelope function methods. 

Although atomistic TB simulations capture the real- 
istic features of a CTAP device and are useful to guide 
CTAP experiments of the future, a typical experiment 
will involve scanning over a large bias range with any 
number of possibilities for straggled donor positions. 
Thousands of large scale TB simulations to understand 
the device will be impractical and time-consuming. For 
this purpose, we have developed an effective 3x3 model 
constructed from a few TB simulations of the system. 
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In the perturbative regime, in which gate voltages cause 
small changes in the A\ state of each impurity, this model 
can provide the same information as the actual TB solu- 
tion of the full atomistic Hamiltonian. The advantage of 
this model is that a large number of voltages and strag- 
gled positions can be explored very rapidly, and correc- 
tion procedures can be determined quickly. By imposing 
symmetry conditions and combining numerical parame- 
terization and analytic solutions, we can even find sets 
of gate biases for various points of the adiabatic path, as 
we demonstrate below. 

Assuming three different donor sites, and a wavefunc- 
tion localized in each donor, we can use a 3 x 3 Hamil- 
tonian describing the system in this 3-state basis. This 
Hamiltonian H e R is of the form, 



El t>LM tLR 

H v \\ = ( £lm* Em tMR 

tLR* tMR* Er 



(3) 



where Ei is the on-site energy of the z-th impurity, and 
Uj is the tunneling matrix element from impurity i to 
impurity j. If the on-site energies of the donors are 
aligned closely, the tunneling elements tLM and tMR can 
be switched alternately from on to off adiabatically to 
effect the electron transfer [39| . 

The idea is to make this model more realistic by cal- 
culating the matrix elements from atomistic TB simu- 
lations. The TB Hamiltonian for a single donor under 
zero gate bias is given by, Hi = Ho + Vi(Ri), where Hq 
is the crystal Hamiltonian, and Vi is the core corrected 
Coulomb potential of impurity i = L,M,R located at 
Ri = Rl,Rm,Rr- If we solve for Hi for each of the 
three donors separately, we can obtain a set of donor 
states {^ji}, where j is the state index, and i is the im- 
purity index. In this case, we restrict our attention to 
the donor ground state only. So, j = 1, and i = L, M, R, 
and the basis set <p = {i/jil, ^im, ^ir}- The full TB 
Hamiltonian of the CTAP device is, 



H T = H + V L (R L ) + Vm(Rm) + V R (R R ) 

+V G (Vsi,V B1 ,V B2 ,Vs2) 



(4) 



The matrix elements of H e R is then given by, [i? e ff]zfc = 
(ipu\HT\ipik)- Whereas the full solution of Ht is time in- 
tensive, requiring many matrix vector multiplications in 
Lanczos, it is fairly fast to evaluate each matrix element 
of the 3x3 model as it involves one matrix vector multi- 
plication between the TB Hamiltonian and a vector, and 
one inner product between two vectors. Since the basis 
set $ is non-orthogonal, the 3x3 problem to be solved 
is of the form, iJ e ff0 = ES(j), where S is the overlap ma- 
trix with elements Sik = (^lzl^ifc)- The eigen solution of 
the i^ e ff gives the molecular states of the CTAP device 
in the basis of the A\ states of the three donors. It is 



to be noted that by simply including the donor excited 
states in the basis (for j > 1), the model can be made 
more comprehensive and applicable to other donor re- 
lated problems as well. It is also fairly easy to extend 
this method to any number of donors. 

Our approach at this point is to find the functional 
dependencies of the matrix elements on the four gate 
biases Vg, and atomic straggle as captured in the donor 
potential Vi(Ri). In our case, we were able to perform 
linear fits of all the elements to the gate biases, with the 
function, 



[H e s]ik = aiVsi + cl 2 Vbi + a 3 V B 2 + a 4 V S 2 + a 5 (5) 

where the coefficients a$ are constants determined from 
the fit. Fig. [8][a) shows some diagonal elements as a 
function of Vbi, while 8(b) shows the same for the off- 
diagonal elements. 
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FIG. 8: (a) Diagonal elements (El, Em, Er) of the 3 x 
3 matrix as a function of Vbi- The other gates are set to 
zero. The donor M is equidistant from L and R. (b) The off- 
diagonal elements (tLM and tMR) with Vbi- (c) The state E2 
obtained in TB from the bias set (0, -0.1128, -0.0903, 0) V. 

Once this parameterized 3 x 3 is obtained, we verify 
that the model is capable of deducing the voltages for the 
adiabatic path. To obtain the midpoint of the path, we 
impose the conditions tLM = ^Mi?, El = Em = Er, and 
Vsi = Vs2 = 0, and obtain the voltages, Vbi = —0.1128 
V and Vbi — —0.0903 V from the matrix elements of 
i^eff. These voltages are offset minutely from our brute 
force simulations by 0.01325 V, 0.0028 V, and 0.0197 V 
for V^i, Vbi, and Vbi, respectively. Performing the ac- 
tual TB simulations with these voltages, we are able to 
obtain the midpoint of the adiabatic path, as shown by 
the wavefunction symmetry in Fig. [8fc). In fact, we 
obtained an improvement over our previous voltages, of 
using only two gates instead of three. Therefore, the 
model is helpful for rapid evaluation of straggle effects. 
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V. CONCLUSION 

In this work, we investigated the effect of imperfec- 
tions on coherent electron transport by adiabatic passage 
in a test triple donor system. Donor misplacements can 
severely hamper the adiabatic transfer path if the mis- 
placements are along the chain, but can be retuned in 
theory, which is a central conclusion of this paper. We 
also performed a study of the sensitivity of the path to 
gate voltage fluctuations showing non-negligible sensitiv- 
ity. We developed a quick and easy way to bridge between 
large scale tight-binding simulations and possible CTAP 
experiments. The effective 3x3 model that we devel- 
oped from TB simulations for this purpose also unites 
the atomic scale precision of TB with intuitive simplicity 
of a toy model. The model in conjunction with tight- 
binding can help guide experiments even in the presence 
of imperfections such as donor straggle and gate voltage 
fluctuations. Furthermore, the model provides a general 
translation of intensive atomic scale calculations to sim- 
plified models used to describe problems in many solid- 
state and quantum-optics applications. 
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